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Abstract 

The transverse momentum distributions of identihed hadrons produced in Pb-Pb cohisions at 
the Large Hadron Cohider (LHC) are studied in the low and intermediate range for px < 5 GeV/c. 
Ah four spectra (vr, K,p, A) can be weh reproduced in the recombination model based on a common 
thermal parton distribution of light and strange quarks and on shower partons emitted in hard 
and semihard jets. Two essential parameters are adjusted to fit the data, one being the inverse 
slope of the thermal distribution and the other revealing the degree of momentum degradation 
in the medium. Various combinations of thermal and shower parton components are calculated, 
showing the dominance of minijets. The effect of minijets is to produce harder baryons than mesons 
resulting in their ratio to peak at around ~ 3 GeV/c. Substantial portion of the jet energy 
is found to be lost to the dense medium before the partons emerge at the surface to undergo 
hadronization by recombination. 
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I. INTRODUCTION 



Recent data from ALICE on Pb-Pb collisions at ^/snn = 2.76 TeV provide the first 
glimpse of the transverse-momentum {px) distributions of identified hadrons produced at 
the Large Hadron Collider (LHC) Although the pt ranges are from to 3 GeV/c only 
for TT and p, and up to 5 GeV/c for and A", which are very low by LHC standards, 
there are already features that are unexpected by extrapolation from RHIC [2] . The proton 
distribution seems to be pushed harder to higher pt so that the p/n^ ratio, as well as the 
p/tt~ ratio, continue to rise up to the highest detected momentum at ~ 3 GeV/c. The 
A/K ratio does show a well-developed peak at p-r = 3 GeV/c, and the magnitude of the 
ratio increases with centrality, exceeding L5 at 0-5% centrality. Since no model-calculation 
has been cited in Refs. l|, |2| to indicate the existence of any satisfactory explanation of the 
data, it is important to investigate this problem and develop a model that can accommodate 
all aspects of the pt distributions for vr, K,p and A production up to 5 GeV/c. 

The region of low and intermediate transverse momenta, < px < 5 GeV/c, is difficult 
for theoretical treatment that can claim validity throughout the whole range. It is too low 
for perturbative QCD and includes pt too high for conventional hydrodynamics. In our 
study here we focus only on the hadronization process for which we use the recombina- 
tion model that has been shown to be effective in interpreting the data at intermediate pt 
jsl. The input involves thermal and shower parton distributions at late time that are phe- 
nomenological in the sense of having some parameters not determined from first principles. 
However, our approach has the advantage of making transparent the common origin of the 
various components that contribute to the spectra of different hadrons produced. Since the 
recombination model has been successful in explaining the large baryon/meson ratio found 
at RHIC j4-6|, it is sensible to extend that approach here in examining its applicability at 
LHC. It is generally recognized that at high energy in the TeV regime minijets are copiously 
produced from semihard scatterings. The density of shower partons from such minijets [3, Isl 
is therefore expected to be much higher than at RHIC. On the other hand, thermal partons 
at hadronization should not depend sensitively on the collision energy, since the hadroniza- 
tion of a system of higher initial density occurs later when the thermal medium can expand 
to a lower local density in order for confinement to take place. Thus the mixture of shower 
and thermal partons at late time at LHC should be very different from that at RHIC. Since 
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shower partons are harder than thermal partons, we therefore expect the produced hadrons 
to be also harder — more so for the baryons than mesons due to their larger number of 
constituents in the recombination model. This is qualitatively the short explanation for why 
the ALICE data exhibit harder baryon spectra compared to the RHIC data. The essence of 
this work is to carry out the investigation quantitatively. 

The dense partonic medium created at LHC can lead to severe attenuating effects on 
semihard partons that traverse the medium. One can focus on the various mechanisms of 
energy loss in QCD, as in [qI, but such studies of medium effects do not lend themselves 
readily to the analysis of the hadronization process outside the medium, which is our goal. 
A template for a phenomenological description of the medium effects has been found in Ref. 



10l | that takes geometrical properties of Au-Au collisions at RHIC into consideration. We 



shall use the general form of that study for LHC, but not its details. With the assumption 
of the dominance of gluon jets for the determination of all shower partons, the effect of jet 
quenching is common to all flavors of shower partons that undergo hadron formation outside 
the medium. Thus the central issues in this work are on the various ways in which thermal 
and shower partons can recombine at pt < 5 GeV/c in Pb-Pb collisions at 2.76 TeV. To 
achieve a satisfactory description of the spectra of all identified species is our immediate 
goal so that the study of hadron production at higher px can follow at a later stage on a 
firm footing. 



II. MESON AND BARYON INCLUSIVE pr DISTRIBUTIONS 



Let us first collect here the basic equations in the recombination model developed previ- 
ously [8|, [lyj . For the invariant p^ distributions of meson and baryons averaged over rj and 
at midrapidity, we may use the one- dimensional form 
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(2) 



where pi is the transverse momentum (with the subscript T omitted) of one of the recom- 
bining quarks qi, whose evolution from a parton emitted at the medium surface to a valon 
(the constituent quark of the hadron to be formed) is a dressing process that preserves its 
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momentum IjJ. The recombination functions (RFs), R'^'^, that include the effects of dress- 
ing and hadronic structure, have been determined previously and will be given below. We 
note here that our dN^/dpT should be identified with the experimental dN/2nd'r]dpT, which 
integrates over all (p, while our quantity is defined as averages over (p and rj. The parton 
distributions can be partitioned into various components represented symbolically by 

^9152 = rr + rs + ss, (3) 
^12,3 = TTT+rrs + rss + sss, (4) 

where T and S are the invariant distributions of thermal and shower partons, respectively. 
There can be many terms within each of the components shown above. 

It should be noted that the nature of time evolution of the system is not specified in 
Eqs. (II]) - dl]). The recombination model describes hadronization at late time. Thus T and 
S distributions are model inputs at that time that are appropriate for particular collision 
systems. The thermal distribution should not depend sensitively on the collision energy, 
since the density of the thermal system must be reduced by expansion to around the same 
level at any collision energy in order for hadronization to occur. Thus the thermal parton 
distribution that we adopt has the same form as before 8j 

dpi 

where C has the dimension of inverse momentum. The prefactor pi is necessary to yield the 
exponential behavior for the thermal component of the hadronic distribution dNh/p^dpT, as 
we shall see after specifying the RF. [See Eq. (l3Tll below.] On the other hand, the properties 
of shower distribution do depend strongly on the collision energy, not only because of the 
increased rate of creation of hard partons, but also due to the quenching effect of the denser 
medium. With ^ used as a parameter that describes an aspect of the momentum degradation 
in the medium, Fi{q,^) is the hard or semi-hard parton distribution at the surface of the 
medium to be discussed in detail below. For now we just show the shower distribution as 

S{P2,0= [ (6) 

where Si{z) is the shower parton distribution (SPD) in a jet of type i with momentum 
fraction z. In Refs. 0, lisl SPD is determined by regarding the fragmentation function (FF), 
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being the recombination product of two shower partons 

xDnx) = I (r^)}^;-(^i'^2,x), (?) 

where the curly brackets denote symmetrization of the leading parton momentum fractions 
xi and X2. 

Using Eqs. ©-([I]) in ([T]), we obtain for pion production 



[ dpidp2e-^P^^P^^/^R^pi,p2,PT), (8) 
dpT J 

I diPii^M I ^ J]F,(g,Or:S(g,pT), (9) 



dpT 
dpT 



J d^P{^,<j>,b) J ^^F,(g,0^(g,PT), (10) 



where 



TS{q,pT) = J ^^nPi)^'(^)^;-(Pi>P2,PT), (11) 

= jDUpT/q). (12) 

P(^, 0, b) is the probability for the dynamical path length to be ^ for a path at angle 
initiated at (xo,?/o), weighted by the nuclear overlap function, and integrated over all 
{xo,yo). Geometrically, depends on where the trajectory is between the center and the 
periphery of the overlap, and dynamically it depends on the energy loss along the path, as 
will be discussed in more detail in the next section. 
For notational brevity we define 

Fi{q,0 = J d^Pi^,<P,b)F,iq,0 (13) 
where ^ depends on and b, and get from Eqs. ([9]) and (fTTj) 

dN^ = ^ / ^^r%p,)s'^{p„0R:M,P2,PT), (14) 

PrdpT Pt J Pi P2 

where 

5^(P2, = JjYl ^^(^' 0S!{p2/q) (15) 



that follows obviously from Eqs. ([6]) and (1130 . except to emphasize that the medium effect is 
in Fi{q, not in S^{p2/q). That is, the FF in Eq. ([7]) is for partons outside the medium and 
is therefore not modified by it. The factor of 2 in Eq. (I14p arises from the two distinguishable 
components of T"(l)5"'(2) + T'^(l)iS"(2). Such a factor is absent in the TT term because 
T"(1)T'^(2) and T'^(1)T"(2) are indistinguishable in a thermalized medium of q and q. For 
K'^ production we have for TS two terms of us, so we obtain 

^ = ^ / [r^(Pi)5^(P2,0 + r^{P2)s^{PuO] R''{puP2,pt), (16) 

PrdpT P^Pt J Pi P2 
where S'^ is as defined in Eq. ( !T5|) . but with replaced by Because of the mass difference 

between the two constituents of K, {pi,P2,Pt) is not symmetric under the interchange 
of pi and p2', moreover, the two terms in the square brackets are different. Since we consider 
only hadron production at mid-rapidity here, we may approximate p^ by the transverse mass 
rriT, where 

mT = {ml+piy/'. (17) 

For h = K,p and A, it can be significantly different from px when pt is small, and can cause 
the low-pT spectra to deviate from the exponential behavior of the pion spectrum. 

The recombination of shower partons from the same jet is equivalent to fragmentation, 
so we have 



— (18) 

tpt J q ^ q \q J 



PxdpT rriTPT 

We shall denote it as {SSyK It is also possible for shower partons from adjacent minijets 
to recombine; we shall refer to such processes as {SSY^ . 

For baryon production the parton distribution shown in Eq. (jl]) has many components 
that can be expressed more explicitly, though still symbolically condensed, as 

^719253 = TTT + rrS + nSS)^' + (SSS)^' + T{SSf' + {{SS)^'Sf' + (SSS)'' (19) 

in obvious notation except for {{SSy^Sy^ which means that two shower-partons are from 
1-jet and the third one from a second jet. We have found that all terms from multi-jets are 
small for pt < 5 GeV/c, so in this paper we ignore 2j and 3j and omit the specification Ij, 
since only 1-jet is considered. If we use the simplified notation to abbreviate Eqs. ( IT^ and 
(US]) as 

{TS)^ = 2T,S,, {TS)k = T,Ss + TsS,, (20) 



then in a similar way we shorten the expression for p and A as 



(TTS), 


— TqTq{Su + Sd), 
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— SuSuSd, 
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{TTS)a 


= Tq(TqSs + 2TsSq), 


(24) 


{T{SS))a 


— '^TqSqSs + TgSuSdj 


(25) 


(^SSS)\ 


— SuSdSg, 


(26) 



The recombination functions (RFs) have been determined from the study of hadronic 
structure j^, |ll|, [l^. They are given below in terms of which is the momentum fraction 
of a valon (which plays the role in the scattering problem as a constituent quark does in the 
bound-state problem) relative to the pt of the produced hadron; i.e., yi = Pi/pr, 



^"(1/1,2/2) = 2/12/25(2/1 + 1/2 -1), (27) 
i?^(l/i, 2/2) = B'\a + 1,6+ l)yt+'y'+'6{y, + 2/2 - 1), (28) 
i?''(2/i,2/2,2/3) = ^/>B(2/i2/2)""''2/f+'5(2/i + 2/2 + 2/3-1), (29) 

Nb = [B{a + l,a + /3 + 2)B{a + 1, /3 + (30) 
and gf^ is a statistical factor. From previous studies it has been determined that a = 1, 6 = 2 



where 



jl4 - a = 1.75, /3 = 1.05 for proton [l3| and a = l,/3 = 2 for A [jj]. For pion a = b 







because its mass is especially low due to its being a Goldstone boson, so tight binding of 
the constituent quarks leads to broad distribution in yi. 

An immediate consequence of the momentum-conserving 5-functions in the RFs is the 
simplification of the hadronic distributions. We illustrate that by writing out explicitly the 
pion distributions. Substituting Eq. (127|) in ([H]) and dH]), we get 

^/V^T (^2 

Pt^Pt d 

^ = -^ dpme-'^^^S^iPT-PuO, (32) 
PrdpT Pt Jo 

where the pion mass in p^ = rriT is neglected. The SS component can be written in terms 
of FF as in Eq. (|T8|l . for which the properties of RF are already used in relating S* to D in 



Eq. ([7]). Similar equations as the two above can be exhibited for other hadrons, as done in 
Appendix A. Here we show the thermal component of proton 

= N'C^^e-^-'^, (33) 

PrdpT 

where = gg^NBB{a + 2, f3 + 2)B{a + 2, a + /3 + 4), so that we can emphasize, by comparing 
Eqs. (!3T!) and (!33l) . that pion and proton have the same exponential factors; however, the 
latter has a prefactor p^/irtT that arises from the kinematics of recombination, causing the 
TT and p spectra to differ at low pt- Note that the inverse slopes T for both spectra are the 
same. 

The objective in this section to express the hadronic spectra in terms of the parton 
distributions is now accomplished. In the following it is then only necessary for us to focus 
on the latter in order to calculate the former. 



III. THERMAL AND SHOWER PARTONS 

As noted at the end of the preceding section the thermal distributions of p and vr have 
the same inverse slope T, inherited from the thermal partons, whose invariant distribution 
is given in Eq. (j5]). This is the characteristic property of recombination that deviates from 
the hydrodynamical description where Tgiope depends on the hadron mass due to flow effect 
There is apparent disparity between the p and vr spectra in the data; nevertheless, 
because of the prefactor p^/mT for proton and of the dominant resonance contribution to 
pion for < 1 GeV/c, it does not mean that it is impossible to have a good description of 
both spectra using one value of T, as will be shown below. 

We now extend the universality of the thermal parton distribution to include the strange 
sector in hadron production at LHC. In the TeV realm of collision energies it is expected 
that the light and strange quarks are fully equilibrated before hadronization. If that is so, 
then the same thermal distribution given in Eq. ([5]) should be valid for u, d and s quarks, 
the consequence of which can readily be checked in our model calculation. We note that the 
issue is not one that concerns hadronic masses, since the system under discussion consists 
only of thermal partons in local equilibrium. More specifically, before the recombination 
mechanism is applied, we assume that all gluons are converted to m, d and s quarks and 



their antiquarks, which form the initial states of specific channels for hadronization 
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Thermal gluons do not fragment in the usual sense. Even fragmentation of gluon jets is 
treated by conversion to q and q before recombination. 

For shower partons let us begin by first summarizing the subject of momentum degrada- 
tion treated in Ref. jlfl]. At RHIC the Pt^4> and centrality (c) dependencies of the nuclear 
modification factor Raa{pt,4>,c) impose stringent constraints on the dynamical process of 
energy loss. The problem is complicated because of both the geometrical and dynamical 
aspects of the description of the parton's traversal through a non-uniform medium. For 
a given point (xo,?/o) in the transverse plane where a hard parton is created in the initial 
system in a collision at impact parameter b, one has to calculate the geometrical path length 
£{xo, yo, 0, b) of a path at angle and then the medium effect along that path. As the sys- 
tem expands, i becomes longer, but the local density becomes lower, so those compensating 
effects on the net energy loss result in a dynamical path length ^ that can be related to £ 
through an undetermined parameter but without a time-dependent transport description. 
Upon averaging over all creation points (xo,?/o), one arrives at a probability distribution on 
^, denoted by P{^, (f), b), that relates the average parton distribution Fi{q, (p, b) with momen- 
tum q at the medium surface to the distribution, with a definite ^ by a weighted 
average 

Fi{q,<P,b) = J d^Pi^,^,b)F,{q,0. (34) 

This is a general relationship with details all contained in P(^, 0, b). The parton distribution 
Fi{q,^) is related in turn to the distribution fi{k) at the point of creation by 

Fi{q,0 = j dkkUk)G{k,q,0, (35) 

where fi{k) is the parton density in the phase space kdk. G{k,q,^) is the momentum 
degradation function from k to q: 

G{k,q,0=qS{q-ke-^). (36) 

In {l^ the RHIC data 16| on Raa{pt,4',c) are used to constrain the properties of 
P(^,0, 6). It is found that P{C,,(f),b) can be expressed in some scaling form involving the 
scaling variable z = ^/^, where ^{(p,b) = J (i^^P(^, 0, 6). For our purpose here, let us not 
repeat the details of that study for Pb-Pb collisions at LHC, especially when we restrict 
our attention in this work to only the data at most central collisions. We circumvent the 



complications by deriving a simple parametric form as follows. From Eqs. (fT3|) . (1MI) - (136|) 
we have for the average parton distribution at the surface 

Fi{q, = J dm^, 0, 6)gV«/.(ge«), (37) 

which we parametrize for centrality < c < 0.05 in the form 

F,{q, /t) = k{qfMk{q)), k{q) = nq. (38) 



If the details of P{C,,4>,b) discussed in [lO| are used, it can be shown that Fi{q,^) in Eq. 
( 137|) for a relevant value of ^ can be reproduced very closely by Eq. ( 138|) for a corresponding 
value of K. Since in either case an unknown parameter is needed to represent the effect of 
momentum degradation, we choose the latter expression that does not rely on the details 
of P(^,0, 6), and treat k, as the key parameter to fit the data of all hadronic spectra. The 
significance of k is clearly the average momentum fraction, k~^, that a parton retains upon 
traversing the medium. 

IV. PARAMETRIZATION 

For the thermal parton distribution we use the form in Eq. (|5]) and put 

C = 23.2 (GeV/c)-^ (39) 

that is the value determined at RHIC [8]. We shall let T be adjustable to fit the LHC 
data at low pt- Since our aim is to reproduce the pt distributions of all identified particles 
{tt,K,p,A), the use of one parameter T for all thermal partons is not only economical, but 
would indeed be remarkable, if achieved. Unlike in hydrodynamical studies that consider 
mass- dependent flow effect, our approach incorporates the effects of minijets on the hadronic 
spectra at low px, allowing the recombination processes to determine the similar or dissimilar 
behaviors of different hadronic species. 

n 

The parton distributions, fi{k), at creation have been determined in Ref. [17j] for Au-Au 
collision at 0.2 TeV and for Pb-Pb collision at 5.5 TeV. The form used is 

fJk) = K- ^— (40) 
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where K = 2.5, A, B and n are tabulated for various quarks and gluon. For ^/snn = 2.76 
TeV we use logarithmic interpolation between the two energies for inA, B and n. More 
specifically, for gluons we use 

A = 6.2 X 10^ GeV-^ B = 0.98 GeV, n = 6.22. (41) 

Other parameters can similarly be obtained. The quark contributions are found to be 
significantly less than the gluon contribution. To make our calculation transparent, we shall 
assume the dominance of gluon jets and let the summation over parton species in all hard 
parton terms be represented by a factor cr times the contributions from the gluon jets only, 
which we shall calculate explicitly using Eqs. ( HOj) and ( 14T|) for gluon creation at A; > 3 
GeV/c. That is, we approximate Yli-^i by crFg and shall use a = 1.2 in the calculation, 
since the quark jets amount to roughly 20% of the gluon jet. 

Since the parton distribution Fi{q, k) is not reliable for g < 3 GeV/c, its contribution at 
low q is cut off by a smooth function which we take to be 

g{q)= + (42) 

With this cut-off factor we may then let the q integrals in all equations, such as in Eqs. (jH]) 
and ( JTSl) . to start at g = 0. Shower parton with momentum pj has a momentum fraction 
X = Pj/q that peaks at small x; consequently, even for g > 3 GeV/c, the density of shower 
partons at pj < 1 GeV/c can become very high at LHC. At low px we know that thermal 
partons are important. How to separate the shower partons from the thermal partons at low 
Pj is not very well defined. The energy lost by semihard partons as they traverse the medium 
can enhance the thermal motion of the soft partons in the vicinities of their trajectories, 
thereby contributing to a component of the thermal partons that are intimately related to 
the soft component of the shower partons. The effective T that we shall determine includes 
enhanced thermal partons due to the energy loss of the semihard partons. Because of that 
effect the distributions of shower partons that are to be determined from the FF according to 
Eq. ([7]) must be modified at low x so that S{pi) will not exceed T{pi). That modification is 
discussed in detail in Appendix B. The phenomeno logical basis that supports our approach 



to the problem is the success in treating the pion and proton spectra at RHIC from a 
common partonic distribution characterized by a universal T ^, 
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Here we apply the 



same idea to the problem at LHC. 
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To summarize, we have two basic parameters to adjust to fit a large collection of data 
from LHC. Those parameters are T and k that characterize the thermal medium and its 
quenching effects on minijets. Due to the copious production of minijets at LHC, we expect 
that the relative magnitudes of the parton distributions T and S will be different from those 
at RHIC. They contribute to hadronic distributions in the intermediate range 1 < pr < 5 
GeV/c that is hard to quantify in other approaches. Thus our phenomenological treatment 
that employs two free parameters is a worthwhile endeavor, provided that all four hadronic 
spectra can be reproduced. 
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FIG. 1: (Color online) Transverse momentum distribution of pion produced in Pb-Pb collision at 
y^J/VTv = 2.76 TeV. Data are from Ref. [l| for centrality 0-5%. The black curve is the sum of the 
components calculated for TT (blue), TS (red), SS (purple) recombination. The dashed line is a 
fit to the data, whose difference from the solid is attributed to resonance decay. 



V. RESULTS 

We calculate the four hadronic spectra using equations given in Appendix A. In all 
those equations parton distributions involving the parameter ^ are to be replaced by the 
distributions in terms of k, as we have discussed in going from Eq. fl37p to (1381) . We have 
varied T and k to achieve the best fits of the data reported by ALICE for central (0-5%) 
Pb-Pb collisions at 2.76 TeV The results are shown in Figs. 1-4. The values of those 
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FIG. 2: (Color online) Transverse momentum distribution of K produced in Pb-Pb collision at 
y^SiViV = 2.76 TeV. Data are from Ref. [l| for centrality 0-5%. The data points in the region 
Pr < 2 GeV/c are for K~^, while those for pr > 2 are obtained from the ratio /K'^ and the A*^ 
spectrum in [1|, [20]. The black curve is the sum of the components calculated for TT (blue), TS 
(red), SS (purple) recombination. 
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FIG. 3: (Color online) Transverse momentiim distribution of proton produced in Pb-Pb collision 
at ^/sNN = 2.76 TeV. Data are from Ref. [l] for centrality 0-5%. The black curve is the sum of 
the components calculated for TTT (blue), TTS (red), TSS (purple), SSS (green) recombination. 
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FIG. 4: (Color online) Transverse momentum distribution of A produced in Pb-Pb collision at 

n 

y^sjviv = 2.76 TeV. Data are from Ref. [IJ] for centrality 0-5%. The lines are color-coded as in Fig. 
3. 



parameters are 



T = 0.38 GeV, 



2.6. 



(43) 



The color code for the various lines representing the different components are as follows: TT 
and TTT in blue, TS and TTS in red, SS and TSS in purple, SSS in green, and the total in 
black. The agreement between the black lines and the data points for K^p and A in Figs. 
2-4 are evidently very good throughout the whole region where data exist. 

In the case of pion the calculated total falls below the data for < 1-5 GeV/c. The 
dashed (black) line is drawn to fit that region and represents the extra contribution from the 
decay of resonances. Such a component cannot at present be calculated in this formalism, 
because the RFs have been derived from the hadronic structures of the lowest bound states. 
Resonances involve orbital excitations for which the RFs have not been investigated. The 
dominance of resonance production at low pt is a known fact even for meson-proton collisions 
at ^/s = 53 GeV 20[ where p and u contributions to the tt spectra exceed 60% 21| . At LHC 
not only do vector and tensor mesons decay into pions, but various baryon resonances can 
also contribute to additional pions. Since resonance contribution is also present at RHIC, 



but not considered in Ref. 



we put aside the region at px < I GeV/c so as to give emphasis 



to the difference between what we can calculate at LHC and RHIC for px > I GeV/c. 
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What is notable in these spectra is the important role that the shower parton 5* plays 
throughout the px range shown. In all of them TS and TTS components are of the same 
magnitudes as TT and TTT, or higher. In fact, for px > 0.5 GeV/c, TS is larger than 
TT, and TTS is larger than TTT for all hadrons produced. It means that at LHC minijets 
are pervasive and their effects dominate the spectra throughout the low and intermediate 
Pt regions. For pt < 5 GeV/c, SS,TSS and SSS components are all neghgible. Thus 
the traditional fragmentation of jets can be ignored. In the recombination approach to 
hadronization TS and TTS components represent the medium effect on semihard partons, 
which lose energy to the medium and then regain some of the momenta back by coalescing 
with enhanced thermal partons to form hadrons. That they are important for pj- > 3 GeV/c 
is expected, as has been found in the study of RHIC physics [8], but to see them as being 
so important for pt > 1 GeV/c is a new revelation at LHC. 
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FIG. 5: (Color online) Shower and thermal parton distributions at LHC and RHIC. 

To help to understand more clearly how hadron production at LHC differs from that at 
RHIC, we show in Fig. 5 the shower and thermal parton distributions at LHC and RHIC. 
S and T are as defined in Eqs. (B6) for gluon jet and in (5), and shown by red and blue 
lines, respectively. Appendix B contains the details about S at LHC, which is more than an 
order of magnitude higher than S for RHIC. On the other hand, T is the same at low pi for 
both energies, but the exponential decrease at intermediate pi is steeper at RHIC (T = 0.32 
GeV) than at LHC (T = 0.38 GeV). As a consequence, T at RHIC is much larger than S 
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FIG. 6: (Color online) Transverse momentum distributions of vr at LHC and RHIC. 

at pi < 2 GeV/c and they cross over at around pi = 3 GeV/c, whereas S and T at LHC 
are approximately the same for pi < 2 GeV/c. These differences have significant effects on 
the pion spectra. In Fig. 6 we compare the pt distributions of pions arising from TS and 
TT recombination at the two energies. Evidently, TT dominates over TS at pt < 3 GeV/c 
at RHIC, but TS at LHC is larger at all pt- It is then clear how minijets are so important 
at LHC that the pion spectrum become much harder than at RHIC, and it is not because 
of enhanced flow. 

As shown in Appendix A, there are many terms in our calculation. To get a good fit 
of all four spectra with the adjustment of essentially only two parameters, T and k, is a 
significant achievement that provides a transparent picture of the relative importance of the 
various thermal and shower partons at < 5 GeV/c. The quality of our fits can be seen 
from Figs. 7 (a) and (b) where the particle ratios p/n and A/K are shown. For the former 
the excellent agreement with data is partially due to the capability of our formalism to 
reproduce the p and vr spectra accurately and separately for pt > 1.5 GeV/c, and partially 
because we have included the resonance contribution to the pions at pr < 1-5 GeV/c that 
we did not calculate. Since the reliable portion of our calculation is in the intermediate pt 
region above 1.5 GeV/c, the peaking of the p/n ratio at ~ 3 GeV/c and the gentle falloff 
above that are our prediction. The A/K ratio in Fig. 7(b) does not show perfect agreement. 



LHC RHIC 
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but it is in linear scale that amplifies the deviations of the calculated spectra from the data, 
which are seen in Figs. 2 and 4 to be well reproduced in log scale. From that more tolerant 
point of view the general agreement of the calculated ratio with data may be regarded as 
remarkably good. The qualitative notion gained from the data that baryons are produced 
with more transverse momentum 2|] is now given a quantitative interpretation that shower 
partons from minijets harden the pt spectra because of the larger phase space opened up 
at higher px when three quarks recombine, one of which being the harder shower parton (as 
can be seen in Fig. 10 in Appendix B). 
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FIG. 7: Particle ratios of (a) p/n and (b) A/K. The data are from Ref. 



m 



for centrality - 



70. 



ALICE also has data on the pt distribution of all charged particles [23|] . If we regard A° 
as representative of then adding our result of TT,K,p and A should come close to all 
charges. In Fig. 8 we show the sum of those four by the black line, which almost saturates 
the data, leaving very little room for baryons that are not included. The plot ends at ~ 5 
GeV/c where our calculation ends, but the data go on to px ~ 15 GeV/c. The mismatch 
will be large for > 5 GeV/c due to aspects of jet physics that we have not taken into 
account in the present study, but will be examined in the future. 

The value of T at 0.38 GeV is slightly larger than 0.32 GeV that was determined for 
thermal partons at RHIC . That is eminently reasonable for the hotter plasma created at 
LHC. Since the thermal partons that participate in hadronization are at late time, T is not 
the temperature of the system at initial equilibration time, which would be much higher. 
After expansion and cooling the value at LHC is, nevertheless, still higher than at RHIC 
because there are many more minijets produced, which lose energy to the medium and raise 
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FIG. 8: (Color online) Charged particle pt distribution in black line that is the sum of the four 



components vr (red), p (purple), K (blue), and A (green). The data are from Ref. 
0-5%. 



231 ] for centrality 



the ambient T. 

The value of k at 2.6 implies that on average roughly 1 — = 60% of the initial parton 
energy is lost to the medium. Thus the momenta of the hard and semihard partons created 
at LHC suffer severe degradation as they traverse the dense medium. The precise value of 
K, may change at higher pj- since the degree of momentum degradation can depend on the 
initial parton momentum. 



VI. CONCLUSION 



We have studied the properties of hadron production in the difficult intermediate pt 
region that is too low for reliable calculation in pQCD, and too high for hydrodynamical 
description. Although we have been successful in reproducing the data, our main point is, 
however, not so much on the accuracy of the fits of the data, but on learning about the 
physics at LHC. Since our treatment of the problem is focused only on hadronization at 
the final stage of the evolution of the partonic system, we cannot address issues related to 
the dynamical development in time. Nevertheless, the investigation provides quantitative 
information about minijets that are crucial in the understanding of the nature of hadronic 
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spectra. A distinctive feature about our approach is that the production of mesons and 
baryons with widely different masses can be described on the basis of a common system 
of partons without any specific reference to flow. In the intermediate pt region between 
2 and 5 GeV/c, the non-fiow part of the hydrodynamical description can presumably be 
related to the thermal-shower component of recombination. Hadron masses do not enter 
explicitly in our calculation, although they are implicitly involved in the determination of 
the recombination functions. The relative magnitudes of the TT, TS and SS components for 
mesons (and of TTT, TTS, TSS, and SSS components for baryons) provide a clear picture 
of the smooth transition throughout that pr region, and unify the very different hadronic 
sectors. That picture offers an illuminating complement to the conventional hydrodynamical 
description without contradicting it. 

While the shower partons from minijets play the important role in hadronization in the 
intermediate pt region, harder jets will, of course, become more important at higher p^- We 
do not expect our description to change significantly when we extend our consideration to 
that region, except for the need to study also multi-jet recombination. Before doing that, 
we still have to examine the centrality dependence in the intermediate region. That subject 
is interesting not only because of the data on the second harmonic, V2{pt), that is usually 
interpreted as elliptic flow, but especially because the dominance of minijets implies the 
importance of non-flow. The dependence of momentum degradation at RHIC has been 
considered in our approach previously [lo|], and can be adapted for collisions at LHC That 
will provide another test of the validity of our finding here about the minijets. 
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Appendix A: Hadronic Distributions for 7r,if,p and A 

We summarize in this Appendix the formulas for all the hadronic distributions. They 
follow from the general equations described in Sec. II and the specific recombination func- 
tions for the four hadrons. We shall leave out the multi-jet (2j and 3j) contributions because 
they are small and negligible for pt < 5 GeV/c, but we exhibit explicitly the various T 
and S combinations. As explained in Sec. IV we approximate the sum over all hard parton 
species by the dominant gluon jet contribution multiplied by a constant factor a (set to be 
(7 = 1.2) that roughly accounts for the combined contribution of all other quark-jet terms. 
That significantly simplifies the calculation without compromising the general features of 
the results. 



1. Pion 



^ = ire-^^/^, (Al) 
dNf 1 r dq 



PrdpT Pt 

where 



J ^aF,(q,OD';(pT,q), (A3) 



S'^iP2, = / ^^F,{q, OS'lip,, q) (A4) 

The variable q is the momentum of gluon jet at the medium surface, while the index q refers 
to light quarks to be distinguished from s quark whose role in the RF of kaon is different. 
The SPD S"i and FF D'^ are modified at low momentum in ways discussed in Appendix B. 



20 



2. Kaon 



PrdpT 
PrdpT 



NKB(a + 2,b + 2)^e-^^/^ 
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K 



NkC 



dpipl{pT-pif 
X [p.e-^^'-'S^pT - PuO + (PT - Pi)e~^^^-^'y^S^ip,, 01 



(A5) 

(A6) 
(A7) 



PrdpT iii'T 

where S^{p2,i) is as defined in Eq. (A4), except for S'g{p2,q) replacing S'g{p2,q)- The RF 
for kaon is found in Ref. [l4| to be for a = 1 and 6 = 2. 



3. Proton 



where 



dN^'^^ 

PrdpT 
dN^ 

PrdpT 



dN^'S 
PrdpT 



dN^'^ 
PrdpT 



gltNpB{a + 2, /3 + 2)B{a + 2, a + /3 + 4) 



rrir, 



VT 



dpi 



VT-Pl 



dp2{piP2)°'^^ {Pt - Pi - P2y 



X 



p 2a+/3+3 

rrirpprp 



dpi 



PT-Pl 



dp2{PlP2T{pT -Pl-P2y 



Uln 



xpie ^'/^5«''(p2,Pr -Pi -P2,6, 
^aF,{q,OD'lipT,q), 



S''{P2.P,.0 = / '^^F,{qri)S'l{p2.q)S'l{p,,q-p2). 



(A8) 

(A9) 

(AlO) 
(All) 



(A12) 



The factor gf^^ is 1/12 when feeddown from A is excluded jj]. The RF parameters for proton 
are a = 1.75 and /3 = 1.05 Q. 
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4. A 



'T 



dNT^ 
PrdpT 



dNps 
PrdpT 



dNf^ 



gftNAB{a + 2,(3 + 2)B{a + 2, a + /3 + 4) 



mU 



(A13) 
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PT 
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PT-Pl 
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xe 



-{Pl+P2)/T\r)cq 



[2S'^{PT -Pi- P2, + S'iPT -Pi- P2, 0], (A14) 



g^tNAC 



A 2a+/3+3 

mifprp 



Pt 



PT-Pl 



dp2{piP2T{PT -P1-P2Y 



dpi 

xp,e-P^/^[2S'^'^{p2,PT-Pi -P2,0+S''iP2,PT-Pi -P2,0], (A15) 

. ^ / ^cTF,{q,OD'^{pT,q), (A16) 

PrdpT J T 

where ^f^^ = 1/8 (1/2 for A° or and 2/8 from spin consideration). The RF parameters for 
A are a = l,/3 = 2 i5^'^(p25P35 defined in Eq. (A12), but with ^"^(ps) replacing 
S'1{P,)- 



Appendix B: Shower Parton Distribution 

We derive in this Appendix the shower parton distribution we use for the study of hadron 
production at LHC The basic idea is aheady described in Refs. [3, [3]) only the parametriza- 
tion is now different. Fragmentation function (FF) at high in deep inelastic scattering 
gives the hadron distribution in a quark or gluon jet. It does not specify the way in which 
hadrons are formed. In the recombination model FF is described as a two-step process, first 
the development of shower partons in a jet, then the coalescence of the shower partons to 
form a hadron. For pion production it is as expressed in Eq. ([7]), which should be augmented 
by the dependence of Z^.f(x,/i^) and S^{x,fi'^) on the energy scale u. Since the evolution 



of D'^{x,jj?) in /i^ can be tracked experimentally and theoretically |24|, the /i^ dependence 
of can be determined accordingly by use of Eq. ([7]). However, to include that de- 

pendence in the application of Si (x) in Eq. (IT^ for heavy-ion collisions is too complicated 
and more meticulous than necessary in view of the many other unavoidable approximations. 
Thus a fixed yU^ is used in practice. For RHIC we have used /i = 10 GeV [T], [12|]. Now for 
LHC we continue to use the same /z, since the hadronization scale at late time is the same. 
However, we improve the determination of the SPD and include scale-breaking effects due 
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to a cut-off at low pi. 
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FIG. 9: Fit of KKP fragmentation function D'^{x,^^ 



2j] (squares) by Eq. (Bl) using tlie recom- 



bination formulas Eqs. (7) and (27), as shown by the solid line. 



The gluon FF D'^{x,^'^) for /i = 10 GeV is shown by the square points in Fig. B-1, 



24| . To reproduce that x dependence, we use 



following the parametrization given in Ref. 
the SPD from gluon to light quark q 

S^g{x) = ax\l - xYil + dx^) (Bl) 

where a = 0.739,6 = -0.28, c = 4.387, ci = 4.502 and e = 10.469. The solid line in Fig. 9 
shows the result of our calculation based on Eqs. and ( IBll) . The fit is evidently very 
good. 

In the application of S^{x) to hadronization processes in heavy-ion collisions, we have 
discussed in Sec. IV how the shower partons at very low momenta pi are merged into the 
region dominated by thermal partons. The semihard partons lose energy to the medium 
whose thermal partons are enhanced in such a way as to make the distinction between 
thermal and shower partons meaningless. The peaking of the SPD at very low pi is therefore 
unrealistic. Since the inverse slope of the thermal parton distribution is to be determined 
phenomenologically that includes the effect of energy loss of the semihard partons, we require 
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that the SPD should not exceed the exponential behavior of the thermal partons. The scale 
invariant thermal distribution T(pi) given in Eq. ([5]) is shown in Fig. 10 by the red line; it 
dips at small pi because of the prefactor Cpi. The integrated shower distribution S^{pi), 
defined as in Eq. f[T^ but for gluon jet {i = g) only and with ^ replaced by n in Fg{q, k) 
given in Eq. (155]) . is shown as a function of pi by the dashed line Fig. 10. We see that 
exhibits a power-law behavior for pi > 1 GeV/c that is expected, and is larger than the 
exponential behavior of T{pi). However, the peaking at low pi for pi < 0.5 GeV/c is due to 
the unreliability of the method of determining S^{x) at low x. Since for physical reasons we 
want the thermal partons to dominate at low pi, we introduce a cut-off factor on the SPD. 
A cut-off with a particular scale implies breaking of scale-invariance described by x. Such 
a breaking is reasonable in low-g processes. If we write x = Pt/q in Dg{x), but now in the 
non-scaling form Dg{pT, q), we introduce the low-pr cut-off as follows 

D'lipT, q) = D-gipT, qhiipT), 7i(Pt) = 1 - e'^'-. (B2) 
The shower parton distribution must therefore also be modified in accordance to 

^D'lipT, q) = - dp^S'lipu q)S'l{pT -Puq-Pi), (B3) 
Q Pt Jo 

where Eq. fl271) is used in ([7]) with xi = pi/q. The modified 5"^(pi, q) is now 

S'l{p,,q) = S^g{p,/qh2{pi), (B4) 

where the corresponding cut-off in pi is 

72(pi) = l-e-(^'i/°-3)'. (B5) 

It can be demonstrated that the two damping factors 7i(pt) and 72(^1) are coordinated to 
satisfy Eq. (IB3p for all px in the range of 1 < < 15 GeV/c. 

The consequence of this cut-off on the integrated shower distribution 

S'lipuK) = I ^F^(g,«:)^'^(pi,g) (B6) 

is shown by the solid (black) line in Fig. 10. Note that it is now just lower than the red line 
of the thermal distribution T{pi) for pi < 0.5 GeV/c, which is the criterion of this cut-off. 
In actual computation of the hadronic spectra this modified distribution S' is used in all 
shower distributions generically expressed as S in Sec. II and Appendix A. 
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FIG. 10: (Color online) Thermal distribution T{pi) is depicted by the blue line for T = 0.38 GeV. 
Shower parton distributions are shown by Sg{pi) (dashed line) and 5'^(pi) (solid red line); the 
latter includes the low-pi cut-off. 
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